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Improved  Soil  Erosion  and 
SedimentTransportin  GSSHA 


by  Charles  W.  Downer,  Fred  L.  Ogden,  Nawa  Pradhan,  Siqing  Liu,  and  Aaron  R.  Byrd 


PURPOSE:  To  describe  the  new  sediment  transport  fonnulation  in  GSSHA  and  demonstrate 
the  improved  capability  in  predicting  storm-total  sediment  runoff  using  a  research-quality  data 
set. 

BACKGROUND:  GSSHA  simulates  overland  soil  erosion  and  outputs  erosion  and  deposition 
for  any  size  class  of  particles  smaller  than  gravel  with  specific  gravities  greater  than  1.0.  The 
model  first  calculates  particle  detachment  by  raindrops  and  surface  runoff,  and  then  calculates 
the  transport  capacity  of  surface  runoff  using  one  of  three  user-selected  transport  equations.  The 
actual  sediment  transport  is  determined  by  comparing  the  amount  of  detached  soil  and  the  trans¬ 
port  capacity  of  surface  runoff.  Depending  on  particle  size,  soil  transported  to  channels  is  treated 
either  as  wash  load  or  bed  load,  with  sizes  less  than  the  user  specified  value  of  sand  treated  as 
wash  load  and  larger  sizes  treated  as  bed  load.  The  original  GSSHA  sediment  transport  formula¬ 
tion  was  hard-coded  to  simulate  only  three  size  classes  of  sediment  (sand,  silt,  and  clay),  each 
with  a  specific  gravity  of  silicate  minerals  (S=2.65).  The  new  GSSHA  erosion  routines  have  been 
generalized  to  allow  simulation  of  arbitrary  size  classes  smaller  than  gravels,  each  with  unique 
specific  gravity. 

Ogden  and  Heilig  (2001)  perfonned  a  detailed  investigation  of  the  CASC2D  erosion  formulation 
(Johnson  1997,  2000)  and  found  that  it  greatly  over-predicts  stonn-total  sediment  yield  for 
events  that  are  considerably  larger  than  the  event  used  to  calibrate  the  model.  They  attributed  this 
to  the  lack  of  detachment  limits  in  the  original  CASC2D  formulation,  which  assumes  that  trans¬ 
port  capacity  of  overland  flow  is  always  satisfied  by  erosion.  Kalin  and  Hantush  (2006)  com¬ 
pared  the  Kineros-2  model  with  the  CASC2D  erosion  formulation  included  in  prior  versions  of 
GSSHA.  Their  study  revealed  that  while  the  hydrologic  components  of  GSSHA  out-perform 
Kineros-2,  the  lack  of  detachment  limits  in  the  CASC2D  erosion  formulation  was  problematic. 
This  re-formulation  of  the  erosion  routines  in  GSSHA  was  performed  partly  to  address  this  issue. 
Kalin  and  Hantush  (2006)  also  suggested  that  the  raster  channel  representation  in  GSSHA  dis¬ 
torted  channel  lengths  and  slopes,  requiring  the  use  of  artificially-low  values  of  Manning’s  n.  To 
address  this  issue  the  channel  representation  in  GSSHA  was  changed  from  raster  to  vector.  The 
effect  of  this  change  on  representativeness  of  the  Manning’s  n  values  used  in  the  enhanced 
GSSHA  formulation  is  examined  as  well. 

SEDIMENT  TRANSPORT  FORMULATION:  To  address  the  many  concerns  identified  with 
the  original  sediment  transport  formulation  developed  in  CASC2D,  the  sediment  erosion  and 
transport  formulation  in  GSSHA  was  redeveloped  to  function  in  the  following  general  manner. 
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1 .  Soils  on  the  overland  flow  plane  are  comprised  of  a  user  specified  number  of  sediments 
with  specified  diameters,  specific  gravities,  and  fractional  distribution  in  each  overland 
flow  cell. 

2.  GSSHA  keeps  sediment  particles  in  three  separate  bins,  parent  soil  material,  suspended 
sediment,  and  deposited  materials. 

3.  During  a  rainfall  event,  sediment  particles  are  detached  by  rainfall  and  overland  flow. 
Deposited  materials  are  eroded  first,  if  available.  If  there  is  no  deposited  layer,  or  the 
amount  of  deposited  materials  is  less  than  the  detachment  forces,  then  parent  soil  mate¬ 
rials  are  eroded.  This  detachment  process  sets  an  upper  limit  on  the  amount  of  sediments 
that  can  be  transported  from  a  given  cell. 

4.  Next,  the  transport  capacity  is  calculated  using  the  user  specified  method.  If  the  transport 
capacity  is  greater  than  the  available  materials,  then  all  of  the  suspended  materials  are 
transported  to  downstream  cells  by  advection,  using  a  two-dimensional  sequential  expli¬ 
cit  method. 

5.  If  the  transport  capacity  cannot  accommodate  all  of  the  suspended  materials,  then  some 
materials  may  be  deposited  according  to  the  calculated  trap  efficiency.  Materials  depo¬ 
sited  in  cells  are  added  to  the  deposition  layer  in  that  cell  and  are  eroded  first  in  subse¬ 
quent  time  intervals. 

6.  Materials  not  deposited  or  transported  remain  in  suspension,  and  may  be  transported  or 
deposited  in  subsequent  time  steps. 

7.  Upon  reaching  a  stream  cell,  suspended,  and  deposited  materials,  in  an  overland  flow  cell 
containing  a  stream  node  are  added  to  the  stream  node. 

8.  Once  in  the  stream  network,  particles  are  attributed  as  either  wash  load  or  bed  load 
depending  on  whether  the  particles’  diameters  are  larger  or  smaller  than  the  user- 
specified  value  of  sand. 

a.  Particles  sand  size  and  larger  are  transported  as  bed  load  using  Yang’s  method. 

b.  Particles  smaller  than  sand  size  are  considered  wash  load  and  are  transported  using 
the  advection-dispersion  equation.  Particles  considered  wash  load  do  not  settle  in 
fluvial  stream  links  or  structures  but  may  settle  in  reservoirs  depending  on  their  cal¬ 
culated  fall  rate  and  reservoir  properties. 

9.  Sediment  transport  is  now  compatible  with  long-tenn  simulations.  If  long-term  simula¬ 
tions  are  specified  by  the  user,  changes  to  the  parent  and  depositional  layer  on  the  over¬ 
land  flow  plane  and  channel  cross  section  in  one  event  are  carried  forward  to  the  next 
event.  If  the  user  chooses,  changes  in  overland  flow  cell  elevation  due  to  erosion  or  depo¬ 
sition  can  also  be  carried  forward  to  the  next  event.  To  save  computational  time,  any  sus¬ 
pended  sediments  remaining  in  suspension  at  the  end  of  a  specified  event  are  deposited 
into  the  overland  flow  cell  containing  the  suspended  sediments  and  sediment  computa¬ 
tions  cease  until  the  next  specified  event. 
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Sediment  Detachment.  Soil  detachment  can  occur  due  to  rainfall  and  overland  flow.  Total 
detachment  is  comprised  of  the  sum  of  rainfall  and  overland  flow  detachment. 

Detachment  by  raindrops.  Detachment  by  raindrops  is  considered  to  be  a  function  of  rainfall 
momentum,  which  is  also  a  function  of  rainfall  intensity.  Raindrop  detachment  also  takes  into 
account  additional  factors  including  surface  water  cushion,  ground  cover,  and  plant  interception 
(Foster,  1982;  Wicks  and  Bathurst,  1996): 

Dr  =  KlCwCGCiMR  (1) 

where: 

Dr  =  detachment  capacity  rate  (kg  m'  -s'  ) 

Ki  =  soil  erodibility  factor  for  detachment  by  raindrop  impact  (J'1) 

Cw  =  water  depth  correction  factor 
Cg  =  canopy  cover  factor 
C,  =  a  cover-management  factor 
Mr  =  moment  squared  for  rainfall  ((kg  m  s'  )  m'  s'  ) 

Table  1  lists  a  set  of  values  of  the  soil  erodibility  factor  for  different  soils.  These  values  were 
calculated  from  experimental  data  by  various  researchers,  as  listed  in  Table  1. 


Table  1. 

Soil  erodibility  factor,  K/,  for  detachment  by  raindrop  impact,  after  Wicks  and 
Bathurst  (1996) 

Data  Source 

Soil  Type 

Clay 

Silty  Clay 

Silty  Clay  Loam 

Silt 

Silt  Loam 

Loam 

Sandy  Loam 

Sand 

Meyer  and  Harmon(1984) 

19.0 

18.2 

16.2 

29.8 

39.8 

28.2 

32.0 

Morgan  (1985) 

30.0 

Bradford  et  al.  (1987a,b) 

73.5 

22.2 

25.7 

37.6 

34.4 

62.4 

Verhaegen  (1987) 

24.7 

23.4 

30.0 

The  relationship  between  the  momentum  squared  and  intensity  of  rainfall  is  highly  nonlinear.  It 
is  expressed  as  a  power  function: 

Mr  =  alp  (2) 

Where  /  is  rainfall  intensity  (mm  h'1),  and  a  and  (i  are  empirical  coefficients  related  to  rainfall 
intensity.  Standard  values  for  these  parameters  are  shown  in  Table  2. 
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Table  2. 

Parameters  for  the  relationship  between  momentum  squared 
and  rainfall  intensity  /,  after  Wicks  and  Bathurst  (1996). 

Rainfall  Intensity  (mm  h'1) 

a 

B 

0</<10 

2.69*1  O'8 

1.6896 

10</<50 

3.75*1  O'8 

1 .5545 

50</<100 

6. 12*1  O'8 

1 .4242 

1 00</<250 

11. 75*1 0'8 

1.2821 

The  calculation  of  water  depth  correction  factor  is  based  on  the  assumption  that  surface  water 
depths  greater  than  a  critical  depth  protect  the  soil  from  raindrop  impact  erosion.  The  expression 
for  the  correction  factor,  Cw,  is: 

r  (U/i\  fexP0-VDJ  if  h>Dm 

c-(Mn  i  if  k<D.  <3> 

where  h  =  water  depth  (m),  and  Dm  =  median  raindrop  diameter  (m). 

The  median  raindrop  diameter  is  determined  from  the  Laws  and  Parsons  (1943)  equation: 

Dm  =  0.00 1247° 182  (4) 

Detachment  by  surface  runoff.  Surface  runoff  detaches  soil  particles  by  exerting  a  shear 
stress  that  breaks  the  bonds  between  particles.  Erosion  in  rills  is  lumped  and  described  as  gross 
rill  erosion.  Within  a  grid  cell,  rill  erosion  and  flow  are  assumed  to  be  uniformly  distributed.  The 
detachment  capacity  rate  by  surface  runoff  has  the  form: 

Dc  =  a(r-rJ(l-GlTc)  (5) 

where: 

2  1 

Dc  =  detachment  capacity  rate  (kg  m'  -s'  ) 
a  and  b  are  empirical  coefficients 
r  =  the  flow  shear  stress  (Pa) 

Tcr  =  the  critical  shear  stress 
G  =  the  sediment  load  (kg  in'  -s'  ) 

2  1 

Tc  =  the  sediment  transport  capacity  of  surface  runoff  (kg  m'  -s'  ) 

In  the  model,  it  is  assumed  that  detachment  is  linearly  proportional  to  the  excess  shear  stress,  and 
b  is  taken  as  1  (Flanagan  and  Nearing,  1995). 
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Sediment  Transport  Capacity  of  Surface  Runoff.  Two  alternative  methods  have  been 
added  to  the  GSSHA  formulation  to  calculate  the  sediment  transport  capacity  of  surface  runoff  in 
addition  to  the  original  Kilinc-Richardson  (1973)  equation. 

The  Kilinc  and  Richardson  (1973)  sediment  transport  equation,  as  modified  by  Julien  (1995), 
and  Ogden  and  Heilig  (2001)  is  used  in  the  current  GSSHA  code.  Sediment  discharge  by  means 
of  overland  flow  is  a  function  of  the  hydraulic  properties  of  the  flow,  the  physical  properties  of 
the  soil,  and  surface  characteristics  as  given  by: 

q ,  =  25500^2  0355'-664  —  (6) 

f  0.15 

where: 

qs  =  sediment  unit  discharge  (ton  m"1  sec'1) 

2  1 

q  =  unit  discharge  (m  sec'  ) 

Sf  =  friction  slope  (-) 


The  factor  K  in  Equation  6  is  a  combined  factor  that  describes  soil  erodibility,  land-use  and  land- 
cover  characteristics.  It  can  be  conceptualized  to  consist  of  multiplicative  terms  similar  to  the 
USLE  soil  erodibility  factor  (0-1),  soil  cropping  factor  (0-1)  and  conservation  factor  (0-1)  in  the 
development  by  Julien  (1995).  The  use  of  one  factor  A  represents  a  departure  from  Julien  (1995), 
who  used  all  three  factors  from  the  Universal  Soil  Loss  Equation  (USLE).  This  departure  is  justi¬ 
fied  by  questions  regarding  the  applicability  of  the  USLE  terms  at  the  time  scale  of  individual 
storm  events. 


The  Engelund-Hansen  (1967)  equation  can  be  used  to  calculate  sediment  transport  for  each  soil 
size  and  the  resulting  total  transport  is  calculated  by  multiplying  the  proportion  of  the  size  in  the 
parent  material  by  the  calculated  rate 


__  0.05 BV2h3,2S3/2 

Ar  - — - 

(s  -!)"£>  y[g 


where: 


Gi  =  the  volumetric  sediment  transport  rate  of  z'-th  size  fraction 
K  =  the  calibration  coefficient  (=  1  for  standard  equation) 

Fj  =  the  proportion  of  z-th  faction  in  the  parent  material  or  deposited  layer 

B  =  the  width  of  flow 

V  =  the  mean  water  velocity 

h  =  the  flow  depth 

S  =  the  water  surface  slope 

s  =  pjp  is  the  specific  gravity  of  z'-th  fraction 

g  =  the  gravitational  acceleration 


(7) 
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l os  =  the  sediment  density 
p  =  the  water  density 
Di  =  the  mean  size  of  /-th  fraction 

The  factor  0.05  in  Equation  7  was  determined  using  empirical  data.  The  suggested  applicability 
of  the  Engelund-Hansen  equation  is  for  sJdi5  /D25  <1.6  ( Dn  is  the  grain  size  for  which  N  percent 
of  sediment  is  finer  by  weight)  and  for  sand-size  sediments  coarser  than  0. 15  mm. 

Several  excess-shear  methods  that  are  similar  to  those  used  in  the  WEPP  (Flanagan  and  Nearing, 
1995)  model  can  also  be  chosen  to  calculate  the  transport  capacity  of  surface  runoff.  The  general 
expression  has  the  form: 

Tc=  a(a>-CQc)h  (8) 

where: 

a  =  the  transport  coefficient 
b  =  the  exponent 

co  =  shear  stress,  steam  power,  or  unit  stream  power 
coc  =  a  critical  shear  stress,  stream  power,  or  unit  stream  power 

A  comparison  of  the  three  sediment  transport  relations  denoted  by  Equations  6,  7,  and  8  was  per¬ 
formed  by  running  the  GSSHA  model  with  each  equation.  Surprisingly,  the  result  of  this  com¬ 
parison  showed  that  for  sediments  with  S=2.65,  there  is  very  little  difference  between  them. 
Therefore,  the  user  is  advised  to  use  the  Kilinc-Richardson  (1973)  method  because  it  has  the 
smallest  number  of  parameters.  However,  for  simulations  involving  sediments  with  specific 
gravities  different  from  2.65,  the  use  of  the  Engelund-Hansen  (1967)  equation  is  required. 

Sediment  Transport  in  Channels  with  Breakpoint  Cross-sections.  The  GSSHA 
model  employs  the  unit  stream  power  method  of  Yang  (1973)  for  routing  sand-size  total-load  in 
stream  channels.  Unit  stream  power  is  defined  as  the  product  of  the  average  flow  velocity  U  and 
the  channel  slope  S.  Soil  erosion  in  channels  is  considered  transport-limited. 

In  the  GSSHA  model,  the  routing  fonnulation  for  sand-size  sediments  works  for  both  trapezoidal 
and  natural  channel  cross-sections.  The  degradation  of  trapezoidal  channels  is  assumed  to  occur 
uniformly  from  the  bottom  of  the  channel.  Degradation  can  continue  and  bed  load  is  transported 
at  the  rated  calculated  with  the  Yang  (1973)  method  until  the  user-specified  maximum  degrada¬ 
tion  is  reached.  Deposition  is  assumed  to  cover  the  entire  bottom  width  of  the  channel,  whether 
in  the  previously  eroded  bottom  or  in  the  trapezoidal  section. 

In  the  case  of  natural  channels  with  breakpoint  cross-section  definitions,  GSSHA  interpolates 
from  the  known  X-Y  pairs  to  get  60  X-Y  pairs  for  each  breakpoint  cross-section.  Channel  geo¬ 
metry  is  updated  after  each  rainfall  event  that  generates  surface  runoff.  GSSHA  tracks  the  high¬ 
est  water  level  in  the  channel  during  a  rain  event  and  uses  it  as  a  point  below  which  channel 
erosion  or  deposition  is  applied  as  a  linear  function  of  depth,  increasing  from  zero  at  the 
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maximum  water  surface  elevation  to  a  maximum  at  the  channel  bottom.  To  take  the  gravitational 
force  into  consideration,  deposition  is  assumed  to  be  proportional  to  the  square  of  water  depth, 
starting  from  the  recorded  maximum  water  surface  elevation.  The  algorithm  implemented  in 
GSSHA  calculates  the  elevation  change  of  all  points  below  the  highest  water  level  so  that  the 
mass  of  sediment  is  conserved. 

Particles  with  sizes  smaller  than  the  user-specified  value  of  sand  are  assumed  to  be  in  suspension 
upon  entering  the  channel,  and  are  transported  as  wash  load  using  the  advection-diffusion 
method  (Downer  and  Byrd,  2007).  Routing  of  suspended  fines  is  a  natural  extension  of  the  expli¬ 
cit  diffusion  wave  channel  routing  method.  Suspended  fine  sediments  are  routed  as 
concentrations. 

Example  Application  on  Goodwin  Creek.  The  USDA-ARS  Goodwin  Creek  Experimental 
Watershed  (GCEW)  data  set  was  used  to  test  the  modified  GSSHA  model.  The  GCEW  is  a  21.1 
km  watershed  near  Batesville,  Mississippi.  The  USDA-ARS  has  collected  rainfall,  runoff,  and 
suspended  sediment  data  in  the  watershed  since  1986  (Figure  1).  The  goals  of  numerical  tests 
included  (1)  testing  of  the  hydrological  component  after  the  new  representation  of  channel  net¬ 
work  was  introduced;  and  (2)  testing  of  the  new  soil  erosion  routines,  especially  their  ability  to 
predict  sediment  yield  under  heavy  rainfall  events,  where  Ogden  and  Heilig  (2001)  found  that 
the  original  CASC2D  erosion  fonnulation  greatly  overestimated  erosion  during  events  that  are 
significantly  larger  than  the  calibration  event. 


Figure  1 .  Location  of  sampling  instruments  at  GCEW  (after  Downer  and  Ogden  2004). 

The  hydrologic  and  sediment  measurements  from  the  USDA-ARS  National  Sedimentation 
laboratory  Goodwin  Creek  Experimental  Watershed  (GCEW),  which  were  used  to  test  the  soil 
erosion  and  sediment  routine  in  CASC2D  by  Ogden  and  Heilig  (2001),  were  used  to  test  the 
sediment  transport  and  other  new  features  affecting  sediment  transport  in  the  GSSHA  model. 
The  new  features  required  re-calibration  and  validation  of  the  model.  The  Goodwin  Creek  data 
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set  from  May  22,  1982  to  August  30,  1982  (Day  of  year  142  to  243)  was  used  to  calibrate  the 
model  and  data  from  1999  were  used  to  verity  the  hydrologic  response  of  the  model  (Downer, 
2008).  Hydrologic  calibration  included  identification  of  a  suitable  value  for  the  Manning  rough¬ 
ness  coefficient  of  the  channel  network.  In  the  calibration  process,  the  channel  roughness  coeffi¬ 
cient  increased  from  0.028  to  0.035,  reflecting  the  shorter  average  length  of  the  links  (Downer, 
2008).  The  value  of  0.035  is  a  more  physically  realistic  value  for  the  GCEW  channels.  The  num¬ 
ber  of  links  in  the  Goodwin  Creek  channel  network  was  increased  from  18  to  29  to  more  accu¬ 
rately  reflect  the  observed  drainage  network. 


Figure  2.  Observed  and  simulated  discharge  from  GCEW  calibration  period  (after  Downer,  2008) 

Calibration  results  are  shown  in  Figure  2.  Mean  Absolute  Error  (MAE)  for  event  peak  discharge 
and  volume  for  this  period  was  21  and  27  percent,  respectively  (Downer,  2008).  Calibration  of 
the  erosion-related  parameters  focused  on  erodibility  factor  K  for  nine  separate  land  use/soil 
complexes,  and  proceeded  using  the  following  approach:. 

1.  The  Kilinc-Richardson  (1973)  equation  was  used  to  predict  the  transport  capacity  of 
overland  flow. 
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2.  Channel  suspended  load  was  simulated  using  advection-dispersion  as  in  the  original 
GSSHA  formulation,  however  in  the  new  non-orthogonal  channel  network. 

3.  The  value  of  the  raindrop  erosivity  coefficient  K1  in  Equation  1  was  assumed  to  be  equal 
to  60  J-l. 

4.  The  soil  erosivity  values  for  each  unique  soil  texture/land-cover  class  was  identified 
using  the  Shuffled  Complex  Evolution  method  (Duan  et  al.  1992). 

5.  The  coefficient  a  used  for  surface  detachment  in  Equation  5  was  2x10-5. 

Sediment  modeling  results.  Figure  3  shows  the  sedographs  computed  by  GSSHA  during  the 
calibration  period.  Sedograph  numbers  in  Figure  3  correspond  to  the  hydrograph  numbers  in 
Figure  2.  The  sedograph  for  the  small  non-numbered  event  between  hydrographs  2  and  3  is 
included  in  Figure  3.  Table  3  lists  the  observed  and  simulated  sediment  runoff  volumes,  and 
error  in  sediment  volume  for  the  events  shown  in  Figure  2.  One  additional  event  is  included  in 
the  analysis  in  Table  3.  The  hydrograph  from  this  small  event  is  not  pictured  (NP)  in  Figure  2. 


Table  3 

Sediment  runoff  performance  of  GSSHA  during  calibration  period.  Note  event 
numbers  correspond  to  sedographs  shown  in  Figure  2. 

Event 

Observed  Peak 
Discharge  (m3  s'1) 

Observed  Discharge 
Volume  (1000  m3) 

Observed  Sediment 
Volume  (m3) 

Simulated  Sediment 
Volume  (m3) 

Error 

(percent) 

1 

23.8 

224.3 

624 

369 

-41 

2 

12.9 

168.6 

500 

184 

-63 

NP 

0.98 

19.6 

3 

3.8 

22 

3 

15.7 

106.6 

240 

181 

-24 

4 

6.1 

62.2 

30 

96 

220 

5 

18.3 

192.0 

256 

269 

5 

6 

9.9 

77.9 

87 

23 

-73 

7 

5.7 

48.8 

35 

28 

-20 

8 

38 

375.9 

640 

640 

0 

Sum 

2415 

1797 

-26 
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Sedograph  Event  1 


Day  of  Year  1982 


Sedograph  Small  Event  Between  2  and  3 


Day  of  Year  1982 


Sedograph  Event  4 


Day  of  Year  1982 


Sedograph  Events  6,7 


226  226.5  227  227.5  228 

Day  of  Year  1982 


Sedograph  Event  8 


Day  of  Year  1 982 


Figure  3.  GSSHA  sediment  discharge  predictions  and  observations  from  Goodwin  Creek. 


10 


ERDC  TN-SWWRP-1 0-3 
August  2010 

As  seen  in  Table  3,  during  the  calibration  period,  the  mean  absolute  error  in  predicted  sediment 
volume  is  52  percent  and  total  error  for  the  simulation  is  -26  percent. 

The  model  was  validated  against  two  large  events  that  occurred  in  1983.  These  events  produced 
peak  discharges  of  71  m  s'  and  148  nr  s'  ,  approximately  2  and  4  times  larger,  respectively, 
than  any  event  in  the  calibration  period.  All  calibrated  parameter  values  remained  the  same.  The 
initial  moisture  was  adjusted  to  better  match  the  discharge  from  the  first  event.  The  model  results 
are  shown  in  Table  4. 


Table  4 

Sediment  runoff  performance  with  GSSHA  during  large  verification  events. 

Event 

Observed  Peak 
Discharge  (m3  s'1) 

Observed  Discharge 
Volume  (1000  m3) 

Observed  Sediment 
Volume  (m3) 

Simulated  Sediment 
Volume  (m3) 

Error 

(percent) 

1 

71 

788 

1216 

1254 

3 

2 

148 

2996 

5066 

2664 

-47 

Comparison  to  prior  results.  One  of  the  many  concerns  expressed  about  the  previous  for¬ 
mulation  of  the  sediment  transport  in  CASCD  by  Ogden  and  Helig  (2001)  was  that  small  events 
were  poorly  simulated.  This  was  thought  to  be  due  to  several  reasons,  including  a  lack  of  con¬ 
straint  on  the  transport  capacity,  errors  in  initial  soil  moisture,  and  errors  in  hydrologic  predic¬ 
tions.  Figure  4  shows  the  percent  error  in  sediment  volume  for  all  calibration  and  verification 
events.  This  simple  analysis  of  the  errors  in  predicting  sediment  volume  and  the  size  of  the  storm 
event,  defined  by  the  peak  discharge,  shows  that  while  there  is  a  general  improvement  in  predic¬ 
tion  error  with  larger  events,  the  model  is  able  to  simulate  events  of  all  sizes,  including  the  small 
events  included  in  the  calibration  period  as  shown  in  Figure  3,  third  panel.  This  is  in  sharp  con¬ 
trast  to  earlier  findings  by  Ogden  and  Helig  (2001)  who  found  that  CASC2D  was  unable  to 
approximately  simulate  the  sediment  volume  from  smaller  events. 
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Figure  4.  Relationship  between  storm  peak  discharge  and  the  error  in  the 
prediction  of  sediment  volume. 


Another  concern  of  Ogden  and  Helig  (2001)  was  propagation  of  errors  -  that  errors  in  simulating 
discharge  would  lead  to  larger  errors  in  predicting  sediment  discharge.  While  this  is  true  in  a 
general  sense,  the  error  in  predicting  discharge  is  less  than  half  of  the  error  in  predicting  sedi¬ 
ment  discharge.  Since  most  errors  in  sediment  volume  fall  around  the  MAE  50  percent  line 
regardless  of  the  accuracy  of  the  discharge  prediction,  there  is  no  general  trend  between  hydro- 
logic  prediction  error  and  sediment  discharge  error  (Figure  5).  Factors  controlling  sediment  dis¬ 
charge  predictions  in  the  model  extend  beyond  the  ability  of  the  model  to  simulate  discharge  in 
the  stream  channel.  One  data  point  that  stands  out  in  both  Figures  4  and  5  is  Event  5.  Although 
both  the  discharge  and  sediment  discharge  was  very  similar  to  those  of  Event  7  of  the  calibration 
period,  the  GSSFIA  model  severely  overestimated  the  sediment  discharge  for  the  event  even 
though  the  prediction  of  discharge  was  just  as  accurate  as  that  for  Event  7.  The  reason  for  the 
anomaly  is  unknown. 


12 


ERDC  TN-SWWRP-1 0-3 
August  2010 


250 


150 


fa 

W 


50 


►  Calibration 
i  Verification 


0 

0.00  10.00  20.00  30.00  40.00  50.00  60.00  70.00  80.00  90.00  100.00 

%  Error  Peak  Discharge 


Figure  5.  Relationship  between  peak  discharge  error  and  sediment  volume  error. 


The  last  concern  of  Ogden  and  Helig  (2001)  addressed  in  this  document  is  that  the  sediment 
transport  routines  in  CASC2D  grossly  overestimated  sediment  discharge  volumes  from  large 
events.  This  was  thought  to  be  caused  by  the  method  always  satisfying  the  transport  capacity 
regardless  of  erosional  supply  of  sediments.  The  verification  events  listed  in  Table  5  indicate  that 
the  model  as  currently  formulated  can  reproduce  the  sediment  runoff  from  very  large  events, 
even  when  no  such  events  are  included  in  the  calibration  period.  And  as  can  be  seen  in  Figure  6, 
the  model  is  capable  of  reproducing  sediment  volumes  from  events  that  range  in  peak  discharge 
and  sediment  volume  more  than  three  orders  of  magnitude  at  the  GCEW. 


□  Simulated 
♦  Observed 


0.1  1  10  100  1000 
Peak  Discharge  (cms) 


Figure  6.  Observed  and  simulated  sediment  discharges  from  calibration  and 
verification  events. 
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Conclusion.  A  new  sediment  transport  formulation  was  developed  for  the  GSSHA  model.  The 
formulation  can  simulate  erosion,  transport,  and  deposition  of  any  number,  size,  and  specific 
gravity  of  particles  in  a  continuous  manner.  This  allows  the  deposition  and  erosion  of  varying 
materials  to  be  tracked  throughout  a  simulation  period  that  can  be  days,  weeks,  months  or  years 
long.  Additional  enhancements  include  detachment  by  raindrop  impact,  flow  detachment  limits, 
and  the  addition  of  the  Engelund-Hansen  (1967)  and  stream  power  soil  transport  equations.  The 
revised  model  shows  good  results  for  predicting  sediment  runoff  volumes  over  an  extended 
3 -month  period.  When  tested  on  events  much  larger  than  any  events  in  the  calibration  event,  the 
accuracy  of  the  model  predictions  did  not  deteriorate.  When  compared  to  the  results  from  pre¬ 
vious  analysis  (Ogden  and  Helig,  2001)  the  model  is  superior  to  the  previous  formulation  in  rela¬ 
tion  to  three  substantial  concerns:  1)  the  ability  to  simulate  small  events,  2)  the  propagation  of 
errors  in  discharge,  and  3)  the  ability  to  simulate  large  events  outside  the  range  of  calibration. 
The  addition  of  detachment  limits  as  well  as  other  improvements  to  the  formulation  represents  a 
clear  improvement  over  previous  versions  of  the  model. 

ADDITIONAL  INFORMATION:  This  technical  note  was  prepared  by  Dr.  Charles  W.  Downer, 
Research  Hydraulic  Engineer,  Coastal  and  Hydraulics  Laboratory,  U.S.  Army  Engineer  Research 
and  Development  Center.  The  study  was  conducted  as  an  activity  of  the  Regional  Sediment 
Management  work  unit  of  the  System-Wide  Water  Resources  Program  (SWWRP).  For  informa¬ 
tion  on  SWWRP,  please  consult  https://swwrp.usace.army.mil/  or  contact  the  Program  Manager, 
Dr.  Steven  L.  Ashby  at  Steven.L.Ashby@usace.army.mil.  This  technical  note  should  be  cited  as 
follows: 

Downer,  C.  W.,  F.  L.  Ogden,  N.  R.  Pradhan,  S.  Liu,  and  A.  R.  Byrd.  2010. 
Improved  soil  erosion  and  sediment  transport  in  GSSHA.  ERDC  TN-SWWRP- 
10-3,  Vicksburg,  MS:  U.S.  Army  Engineer  Research  and  Development  Center 

h  ttps  ://swwrp.  us  ace.  army,  m  i  1/ 
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